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Individuals are always limited by some inelastic resources, such as time and energy, which restrict them 
to dedicate to social interaction and limit their contact capacity. Contact capacity plays an important role in 
dynamics of social contagions, which so far has eluded theoretical analysis. In this paper, we first propose a non- 
Markovian model to understand the effects of contact capacity on social contagions, in which each individual 
can only contact and transmit the information to a finite number of neighbors. We then develop a heterogeneous 
edge-based compartmental theory for this model, and a remarkable agreement with simulations is obtained. 
Through theory and simulations, we find that enlarging the contact capacity makes the network more fragile 
to behavior spreading. Interestingly, we find that both the continuous and discontinuous dependence of the 
final adoption size on the information transmission probability can arise. And there is a crossover phenomenon 
between the two types of dependence. More specifically, the crossover phenomenon can be induced by enlarging 
the contact capacity only when the degree exponent is above a critical degree exponent, while the the final 
behavior adoption size always grows continuously for any contact capacity when degree exponent is below the 
critical degree exponent. 

PACS numbers: 89.75.Hc, 89.75.-k, 87.23.Ge 


Recent empirical studies of Facebook communication 
networks, scientific cooperation networks and sexual con¬ 
tact networks suggest that individuals’ activities are lim¬ 
ited by the time, funds, energy and other inelastic re¬ 
sources. Thus, individuals exhibit limited contact capacity 
(i.e., individuals can only communicate or interact with a 
finite number of neighbors during a short time) in the dy¬ 
namics of epidemic and behavior. Previous studies have 
proven that limited contact capacity enlarges the epidemic 
outbreak threshold and makes the theoretical prediction 
deviate from simulation results more easily. Unfortu¬ 
nately, a systematical investigation the effects of contact 
capacity on the dynamics of social contagions is still lack¬ 
ing. To fill this gap, we first propose a non-Markovian 
behavior spreading model with limited contact capacity. 
Then, we develop a novel heterogeneous edge-based com¬ 
partmental theory for this proposed model, and verify the 
effectiveness of this suggested theory based on large num¬ 
ber of simulations. Through theory and simulations, we 
find that the final behavior adoption size increases with 
the contact capacity. Strikingly, we uncover a crossover 
phenomenon in which the dependence of the final adop¬ 
tion size on the information transmission probability can 
change from being continuous to being discontinuous. We 
find a critical degree exponent above which the crossover 
phenomenon can be induced by enlarging the contact ca¬ 
pacity. However, the final adoption size always grows con¬ 
tinuously for any contact capacity when degree exponent 
is below this critical degree exponent. Our results help us 
to have a deeper understanding of the effects of contact 
capacity on social contagions, and the developed theory 
could be applied to other analogous dynamical processes 
(e.g., information diffusion and cascading failure). 


* tangminghuang521@hotmail.com 


I. INTRODUCTION 

Humans are the basic constituents of the society, and every 
individual can interact with his/her family, friend and peers. 
These interactions among individuals can induce some inter¬ 
esting collective behavior, such as, spontaneous formation of a 
common language or culture, emergence of consensus about a 
specific issue, and the adoptions of innovation, healthy or mi¬ 
crofinance behavior. Understanding the mechanisms or regu¬ 
larities behind these collective behavior has led to a booming 
subfield of research in complex network science - social con¬ 
tagions, which has attracted much attention in recent years Ch 

i. 

Statistical physics approaches were widely used to investi¬ 
gate social contagions. On the one hand, scientists used these 
methods to analyse large databases of social cont^ions, and 
revealed that reinforcement effect widely exists |0]. The re¬ 
inforcement effect means that individual adopting a behav¬ 
ior is based on the memory of the cumulative behavioral in¬ 
formation that he/she received from his/her neighbors. Cen- 
tola established the artificially structured online communities 
to study health behavior spreading, and found that the rein¬ 
forcement effect ^nificantly increases the adoption of a new 
health behavior ||^ 0- The reinforcement effect also exists 
in the adoptions of Facebook ||3l and Skype JH] services. On 
the other hand, researchers proposed some novel models with 
reinforcement effect to describe the dynamics of social conta¬ 
gions. Among these models, linear threshold model ll9l- [Tll] is 
a famous one, and it is a deterministic model (i.e., a trivial case 
of Markovian process) once the network topology and initial 
seeds are fixed. In this model, an individual will adopt the be¬ 
havior once the current fraction of his/her adopted neighbours 
is larger than a static threshold. The linear threshold model in¬ 
duces that the final behavior adoption size first grows continu¬ 
ously and then decreases discontinuously with the increasing 
of mean degree for vanishing small fraction of seeds. Another 
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more realistic way to incorporate the reinforcement effect is 
whether an individual adopts the behavior should take his/her 
cumulati ve p ieces of behavioral information into considera¬ 
tion IO - IIct . In this case, the dynamics is a non-Markovian 
process, which makes it more difficult to develop an accu¬ 
rate theory. Wang et al proposed a non-Markovian behavior 
spreading model, and found that the dependence of final be¬ 
havior adoption size on information transmission probability 
can change from being discontinuous to being continuous un¬ 
der dynamical or structural parameters perturbation S. 

Recently, scholars found that individuals exhibit limited 
contact capacity (i.e., individuals can only communicate or 
interact with a finite number of neighbors during a short time) 
since the inelastic resources (e.g., time, funds, and energy) 
restrict them to dedicate to social interaction from empiri¬ 
cal analysis lEMl. In Facebook communication networks, 
Golder et al revealed that users only communicate with a 
small number of people even though they have many declared 
friends ll20l] . In scientific cooperation networks, a scientist ex¬ 
changes knowledge with only a fraction of his/her cooperators 
in a paper ll^l23] . In sexual contact networks, individuals can 
not have sexual intercourse with his/her all sexual partners in 
a very short time due to the limitation of morality and physiol¬ 
ogy Researchers have studied the effects of contact 

capacity on some Markovian dynamics (i.e., epidemic spread¬ 
ing) 11251 - 1^ . They found that the epidemic outbreak thresh¬ 
old increases when the contact capacity is limited Iz^ . Mean¬ 
while, each connection (edge) has distinct effective spreading 
probability (to be defined in Sec. Hill), which makes the theo¬ 
retical prediction deviate from simulation results more easily, 
especially in the case of strong structural heterogeneity. 

For the dynamics of social contagions, whether an individ¬ 
ual adopts a behavior behavior or not is determined by the 
cumulative pieces of behavioral information that he/she has 
received from neighbors SS. Once the contact capacity is 
limited, the behavioral information transmission will be lim¬ 
ited, thus further affects the dynamics of social contagions. 
However, a systematic study to understand the effects of con¬ 
tact capacity on dynamics of social contagions is still lack¬ 
ing. In this paper, we try to address how the contact capacity 
affects the behavior spreading dynamics. We first propose a 
non-Markovian behavior spreading model with limited con¬ 
tact capacity, in which each adopted individual tries to trans¬ 
mit the behavioral information to a finite number of his/her 
neighbors. In order to understand, quantitatively, the effects of 
contact capacity on behavior spreading, we develop a hetero¬ 
geneous edge-based compartmental theory. We find that the 
final behavior adoption size increases with the contact capac¬ 
ity. More interestingly, the crossover phenomenon in which 
the dependence of the final adoption size on the information 
transmission probability can change from being continuous to 
being discontinuous. By enlarging the contact capacity, the 
crossover phenomenon can be induced only when the degree 
exponent is above a critical critical degree exponent. How¬ 
ever, the final adoption size always grows continuously for 
any contact capacity when degree exponent is below the crit¬ 
ical degree exponent. The theoretical results from the sug¬ 
gested method can accurately predict the above results. 


The paper is organized as follows. In Sec. |II] we describe 
the behavior spreading model with limited contact capac¬ 
ity. We develop the heterogeneous edge-based compartmental 
theory in Sec. |III] In Sec. IIVI we verify the effectiveness of 
the theory through large number of simulations. Finally, we 
present conclusions and discussions in Sec.lVl 


II. BEHAVIOR SPREADING MODEL 


We consider the behavior spreading on uncorrelated con¬ 
figuration networks ll^ |29t] with N individuals (nodes) and 
degree distribution P{k). We use a generalized model SAR 
(susceptible-adopted-recovered) model ifisll to describe be¬ 
havior spreading on networks. At each time step, each indi¬ 
vidual can be in one of the three different states: susceptible, 
adopted, or recovered. In the susceptible state, an individual 
does not adopt the behavior. In the adopted state, an individ¬ 
ual adopts the behavior and tries to transmit the information 
to his/her selected neighbors. In the recovered state, an indi¬ 
vidual loses interest in the behavior and will not transmit the 
information further. Each individual holds a static adoption 
threshold k, which reflects the criterion (wills) of an individ¬ 
ual to adopt the behavior. 

Initially, a fraction of po individuals (nodes) are randomly 
selected to be in the adopted state (seeds), while other in¬ 
dividuals are in the susceptible states. All susceptible indi¬ 
viduals do not know any information about this behavior, in 
other words, the cumulative pieces of information is zero ini¬ 
tially for all susceptible individuals. At each time step, each 
adopted individual v with k' neighbors randomly chooses 
f{k') number of neighbors due to the contact capacity is lim¬ 
ited, and tries to transmit the information to each selected 
neighbors with probability A. Note that the function f{k') 
represents the contact capacity of v, the larger value of f{k'), 
the more neighbors can receive the information from him/her. 
If/(fc') < k', the contact capacity of individual v is limited. 
Once the contact capacity of v is larger than his/her degree, 
we let he/she transmit information to his/her all neighbors. If 

V transmits the information to u successfully, the cumulative 
pieces of information m that u ever received will increase by 
1, and the information can not be transmitted between u and 

V in the following spreading process (i.e., redundant informa¬ 
tion transmission on the edge is forbidden). If m is larger than 
the adoption threshold k, individual u becomes adopted in the 
next time step. From the mentioned procedures of suscepti¬ 
ble individuals becoming adopted, we learn that the dynamics 
of social contagion is a non-Markovian process. The adopted 
individuals then lose interest in the behavior and enters into re¬ 
covered with probability 7. Individuals in the recovered state 
do not take part in the spreading process. The dynamics ter¬ 
minates once all adopted individuals become recovered. 
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III. HETEROGENEOUS EDGE-BASED 
COMPARTMENTAL THEORY 


Taking all possible values of k into consideration, we can get 
the fraction (density) of susceptible individuals at time t 


The non-Markovian behavior spreading model with limited 
contact capacity described in Sec. |II] makes theoretical pre¬ 
diction from the classical theory (e.g., heterogeneous mean- 
field theory) deviate from simulation results easily. On the 
one hand, in this proposed model, whether a susceptible indi¬ 
vidual adopts the behavior or not is dependent on the cumula¬ 
tive pieces of information he/she ever received. In this case, 
the memory effect of non-Markovian process is induced. On 
the other hand, the heterogeneity of effective spreading prob¬ 
ability for edges increases with the heterogeneity of degree 
distribution, and further enhances the difficulty in developing 
an accurate theory. The effective spreading probability of an 
edge includes two aspects: (1) an edge is randomly selected 
with probability f{k')/k', where k' is the degree of adopted 
individual u; (2) the information is transmitted through the se¬ 
lected edge with probability A. Thus, the effective spreading 
probability of an edge for individual v is \f{k')/k'. 

To describe this process, we develop a novel theory - het¬ 
erogeneous edge-based compartmental theory, which is in¬ 
spired by Refs. ll30l - l^ . The theory is based on the assump¬ 
tion that behavior spreads on uncorrelated, large sparse net¬ 
works. We denote S{t), Ait) and Rit) as the density of in¬ 
dividuals in the susceptible, adopted and recovered states at 
time t, respectively. 

Denoting Ok' {t) as the probability that an individual v with 
degree k' has not transmitted the information to individual u 
along a randomly selected edge up to time t. For simplicity, 
we assume that individuals with identical degrees are the same 
in statistics. In the spirit of the cavity theory, we let individ¬ 
ual u in the cavity state (i.e., individual u can not transmit 
information to his/her neighbors but can receive information 
from his/her neighbors). Considering all possible degrees of 
individual v, the average probability that individual u has not 
received the information from his/her neighbors by time t 


5(t)=^P(A:)s(fc,/). (4) 

k 

Similarly, we can get the fraction of individuals who have re¬ 
ceived m pieces of information at time t 

$(m,/) = ^ P(fc)^(fc, m,/). (5) 

fc =0 

According to the definition of Ok' (/), one can further divide 
it as 

= (k'it) + (6) 

The value of (/), (/), and (/) represents that the prob¬ 

ability of individual v with degree k' is susceptible, adopted, 
and recovered and has not transmitted information to its 
neighbors (e.g., individual u), respectively. 

An initial susceptible neighbor individual u of u can only 
get the information from the other k' — 1 neighbors, since in¬ 
dividual u is in the cavity state. Similar to Eq. (|2l), one can get 
the probability that v has m cumulate pieces of information 
by time t 

T{k',m,t) = (1 - ^yO{t)f-^-^[l - 0{t)r. 

(7) 

We further get the probability of individual v in the susceptible 

K — 1 

^k'{t) ='^T{k',m,t). ( 8 ) 

m —0 

If the adopted neighbor individual v with degree k' trans¬ 
mits the information via an edge, this edge will not meet the 
definition of Ok'{t). The conditions of individual v transmits 
information to u are: (1) the edge connecting them is selected 
with probability f{k')/k' and (2) the information is transmit¬ 
ted through this edge with probability A. Thus, the evolution 
of Ok' (t) is 


where k'P{k')/{k) represents the probability that an edge 
from u connects to v with degree k' in uncorrelated network, 
and {k) is the mean degree. It is straightforward to get the 
probability that individual u with k neighbors has m cumu¬ 
late pieces of information by time t 

^(k,mA) = (1 - Po) Q [0(i)]'=-™[l - 0{t)r. (2) 

The formula 1 — po represents that only individuals in the sus¬ 
ceptible state initially can get the information. From Sec. [Ill 
we know that only when it’s cumulate pieces of information 
are less than k, he/she can be susceptible at time t. Thus, in¬ 
dividual u is susceptible by time t with probability 

K—1 

s{k,t) = '^ (f){k,m,t)- (3) 

m—0 


dOk'jt) Xf{k') ,Au. 
dt k' ’ 


(9) 


If f{k') is larger than k', we restrict that v transmits the infor¬ 
mation to his/her all neighbors [i.e., f{k') = k']. 

According to information spreading process described in 
Sec.in the growth of should simultaneously satisfy: (1) 
the adopted individual v does not transmit the information to 
u through the edge between them and (2) v moves into re¬ 
covered state with probability 7. For the first condition, there 
are two possible cases: the edge between u and v is selected 
with probability f{k')/k' and the information is not transmit¬ 
ted through it with probability 1 — A; the edge between u and 
V is not selected with probability 1 — f{k') / k'. From the anal¬ 
yses above, the evolution of is 


dm) 

dt 




A/(fc') 

k' 


( 10 ) 
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Now, combining Eqs. (l9])- (fT0l) and the initial situations [i.e., 
^*fc'( 0 ) = 1 and = 0 ], we obtain the expression of 

(t) in terms of Ok' (t) as 

4«(f)=7[l-0fe'(f)][^^-l]. (11) 

Utilizing Eqs. (| 6 ]l, (O, Q and (fTTI) . we obtain that 


d0k'{t) 

dt 


A/(fc') 

k' 


K— 1 

[0k'{t) - ^ r(fc',m,f)] 

m—0 




( 12 ) 


According to the model described in Sec. [Ill the densities 
of individuals in adopted and recovered individuals evolve as 


dA{t) dS{t) 

and 


(13) 


dR{t) 

dt 




(14) 


respectively. Eqs. (HI and (fT3]) - (fT4l) give us a complete de¬ 
scription of the social contagions with limited contact capac¬ 
ity. The evolution of each type of density versus time can be 
obtained. 

The densities of susceptible, adopted and recovered indi¬ 
viduals do not change when t ^ oo. We denote R{oo) as the 
final behavior adoption size. To obtain the value of R{oo), 
one can first solve 9k' (oo) from Eq. (fT^ . that is 


K — 1 

9k'{oo) = ^ r(fc',m,oo)] -|- 7[1 - 0fe.(oo)][ 

m—0 


k' 

W) 


-!]■ 

(15) 


Iterating Eq. ([I5]l to obtain 9k'{oo). Then, inserting 9k'{oo) 
into Eqs. (Q)-® to get the values of S{oo) and R{oo) = 1 — 
S{oo). 

Another important aspect we mainly focus on is the con¬ 
dition under which the global behavior adoption occurs. The 
global behavior adoption means that a finite fraction of indi¬ 
viduals adopted the behavior, and the corresponding local be¬ 
havior adoption represents that only a vanishingly small frac¬ 
tion of individuals adopted the behavior. Similar to biological 
contagions, we define a critical transmission probability Ac. 
When A < Ac, the behavior can not be adopted by a finite 
fraction of individuals; when A > Ac, the global behavior 
adoption occurs. Now, we discuss Ac for several different val¬ 
ues of po and K. 

Eor po 0 (i-e., only a vanishing small fraction of seeds) 
and K = 1, 9k'{oo) = 1 is the trivial solution of Eq. ( ffSl) . 
If we change the values of other dynamical parameters, such 
as information transmission probability A, a global behav¬ 
ior adoption may occur. The global behavior adoption oc¬ 
curs only when a nontrivial solution of Eq. ([I5]l emerges [i.e., 
9k' (oo) < 1]. Note that the corresponding fraction of 9k' (oo) 



FIG. 1. (Color online) On strong heterogeneous networks, the final 
adoption size R{oo) as a function of information transmission proba¬ 
bility A for (a) different adoption threshold k and (b) different contact 
capacities c. In (a), black circles (n = 1), red squares (n = 2) and 
blue up triangles (k = 3) are the simulation results for c = 1. In 
(b), black circles (c = 1 ), red squares (c = 2 ) and blue up trian¬ 
gles (c = 3) are the simulation results for k = 2. In figure (a) and 
(b), the lines are the theoretical predictions from Eqs. ® and i fTSl l- 
d. We set other parameters as v = 2.1, 7 = 0.1, and po = 0.1, 
respectively. 


should be taken into consideration. Linearizing Eq. (fTSl) at 
9k'{ 00 ) = 1 and summing all possible values of k' one 
can get the critical information transmission probability 


l{k)G{k) 

(P)_(2_^)(fc)’ 


(16) 


where 


G{k) = 

k' 


k'"P{k') 

WW)' 


Note that Ac is tightly correlated with the network topol¬ 
ogy [i.e., degree distribution P{k)] and dynamical parame¬ 
ters [i.e., contact capacity f{k') and recover probability 7 ]. 
Hubs in heterogeneous networks adopt the behavior with large 
probability. Thus, the value of Ac decreases with degree het¬ 
erogeneity. The value of Ac increases with f{k'), in other 
words, increasing the contact capacity of individuals makes 
the network more fragile to the behavior spreading. We should 
emphasize that 7 also affects Ac /7 (i.e., the effective criti¬ 
cal information transmission probability), which has been ne¬ 
glected in previous studies. The value of Ac /7 increases with 
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FIG. 2. (Color online) The final fraction of individuals in the subcrit- 
ical state ‘I>(ft — 1 , oo) versus information transmission probability 
A for K = 2, c = 1 (black circles) and k = 2, c = 2 (red squares). 
The lines are the theoretical predictions from Eqs. © and (Ull-Gl. 
Other parameters set to be i/ = 2.1, 7 = 0.1, and po = 0.1, respec¬ 
tively. 


7 . If f{k') > k' for every value of k' (i.e., adopted individual 
transmits the information to his/her all neighbors), Eq. ( fTbl) is 
the epidemic outbreak threshold If every adopted in¬ 

dividual only transmits information to his/her c neighbors, we 
can get the critical transmission probability Ac of the model in 
Ref. m. 

For Po — 1 ' 0 and k > 1, we find that 6k'{ 00 ) = 1 is the 
solution of Eq. (fTSl) . However, the left and right hands of 
Eq. (fTSl l can not be tangent to each other at 6k' = 1, which in¬ 
dicates that vanishii^y small seeds can not trigger the global 
behavior adoption lH^ . With the increase of po, different de¬ 
pendence of R{oo) on A occurs for different k. That is, the 
growth pattern of R{oo) versus A can be continuous or dis¬ 
continuous. Through bifurcation analysis ofEq. ( ffSl l. we 
find that R{oo) grows continuously for k = 1, while a discon¬ 
tinuous growth may be induced for k > 1. 


IV. SIMULATION RESULTS 

In this section, we verify the effectiveness of the heteroge¬ 
neous edge-based compartmental theory developed in Sec. Hill 
by lots of simulations. For each network, we perform at least 
2 X 10^ times for a dynamic process and measure the final 
fraction of individuals in the recovered [i?(c»)] and subcriti- 
cal state [<!>(« — 1, cxa)]. These results are then averaged over 
100 network realizations. 

To built the network topology, we use the uncorrelated con¬ 
figuration model according to the given degree distri¬ 
bution P{k) ~ k~'^ with maximal degree kmax ~ v/ZV. 
There is no degree-degree correlations when N is very large. 
The heterogeneity of network increases with the decrease of 
i/. For the sake of investigating the effects of heterogeneous 
structural properties on the social contagions directly, the net¬ 
work sizes and mean degree are set to be = 10 , 000 and 
(k) = 10, respectively. All individuals with different degrees 
have the same contact capacity f{k) = c and recover proba¬ 
bility 7 = 0 . 1 . 

We first study the effects of the adoption threshold k and 



FIG. 3. (Color online) On weak heterogeneous networks, the final 
adoption size R{po) versus information transmission probability A 
for (a) different adoption threshold k and (b) different contact capac¬ 
ities c. In (a), black circles {k = 1), red squares {k = 2) and blue up 
triangles (k = 3) are the simulation results for c = 1. In (b), black 
circles (c = 1), red squares (c = 2), blue up triangles (c = 4) and 
green diamond (c = 8 ) are the simulation results for k = 3. And the 
lines are the theoretical predictions, which are solved from Eqs. (0 
and ( I13l l- U4t . We set other parameters as zz = 4.0, 7 = 0.1, and 
Po = 0.1. 


contact capacity c on the final behavior adoption size i?(oo) 
for strong heterogeneous networks in Fig. [T] We find that 
R{oo) decreases with the increase of k, since individuals 
adopting the behavior need to expose more information. Once 
the contact capacity increases (i.e., c increases), individuals in 
adopted state will have more chances to transmit the infor¬ 
mation to susceptible individuals, thus, the values of R{oo) 
increases. Obviously, the theoretical predictions from hetero¬ 
geneous edge-based compartmental theory agree well with the 
simulation results. 

Another important issue we concern is the dependence of 
R{oo) on A. As shown in Fig. [U for strong heterogeneous 
networks the dependence of R{oo) on A is continuous for any 
values of k and c, and we verify this claim by the bifurcation 
analysis of Eq. (fTSl) . We can also understand this phenomenon 
by discussing the fraction of individuals in the subcritical state 
from an intuitive perspective (see Fig. la. An individual in 
the subcritical state means that he/she is in the susceptible 
state, and the m cumulative pieces of information is just one 
smaller than his/her adoption threshold k. From Ref. IIT^ . we 
know that a discontinuous dependence of R{oo) on A will oc¬ 
cur only when a large number of those subcritical individuals 
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FIG. 4. (Color online) The final fraction of individuals in the subcrit- 
ical state $(k — 1 , oo) versus information transmission probability 
A for K = 3, c = 1 (black circles) and ac = 3, c = 8 (red squares). 
The lines are the theoretical predictions from Eqs. © and (Eli-Gl. 
Other parameters are v = 4.0, 7 = 0.1, and po = 0.1, respectively. 



adopt the behavior simultaneously at some information trans¬ 
mission probability. Fig.|2]shows the final fraction of individ¬ 
uals in the subcritical state <i)(K— 1, 00 ) versus A. We find that 
<h(K — 1 , 00 ) first increases and then decreases gradually with 
A, since the existence of strong degree heterogeneity makes 
individuals in the subcritical state adopt the behavior consec¬ 
utively. In these cases, a continuous growth of R{oo) versus 
A occurs on strong heterogeneous networks. 

We now study behavior spreading on weak heterogeneous 
networks, such as i/ = 4.0 in Fig. [3 Similar with the case of 
1 / = 2.1, increasing k leads to the decrease of i?(oo); and the 
value of R{oo) increases with c, that is the network will be¬ 
come more fragile to the behavior spreading once the contact 
capacity increases. Once again, our theory can predict the so¬ 
cial dynamics very well. For the dependence of i?(oo) on A, 
a crossover phenomenon transition is observed. A crossover 
phenomenon means that the dependence of R{qo) on A can 
change from being continuous to being discontinuous. More 
specifically, as shown in Fig. [3b), the dependence of R{oo) 
on A is continuous for small values of c (e.g., c = 1), while 
the dependence is discontinuous for larger c (e.g., c = 8). 
We justify this claim by the bifurcation analysis of Eq. ([I5]l 
from the theoretical view, which is also verified through an¬ 
alyzing (1 )(k — 1, 00 ) from an intuitive perspective in Fig. |4| 
For weak heterogeneous networks, most individuals adopt the 
behavior with the same probability since they have similar de¬ 
grees. When c = 1, $(«; — 1, 00 ) increases continuously with 
A, which leads to a continuous growth in the value of R{oo). 
When c = 8, 4>(ac — 1, 00 ) first increases with A, and reach a 
maximum at some values Ac, and a slight increment of A in¬ 
duces a finite fraction of 4>(«;, c») to adopt the behavior simul¬ 
taneously, which leads to a discontinuous jump in the value of 
R{qo). 

We further study the effects of z/ and A in Fig. |3 for dif¬ 
ferent values of c. For small contact capacity [i.e., c = 1 in 
Figs. 13 a) and (b)], the dependence of R{oo) on A is always 
continuous for any value of In other words, this dependence 
is irrelevant to the network topology. For large contact capac¬ 
ity [i.e., c = 8 in Figs.|3c) and (d)], there is a crossover phe¬ 


FIG. 5. (Color online) The final behavior adoption size versus in¬ 
formation transmission probability and degree exponent, (a) and (c) 
represent, respectively, the color-coded values of R{oo) from nu¬ 
merical simulations in the parameter plane X-u for c = 1, k = 2 and 
c = 8 , At = 3. The theoretical predictions for c = 1, k = 2 and 
c = 8 , K = 3 are shown in (b) and (d), respectively. And theoret¬ 
ical predictions are solved from Eqs. and jl3b - ( I14l l. In (b) and 
(d), R{oo) grows continuously with A to a large value in region I. 
And in region II, R{oo) grows discontinuously and a finite fraction 
of individuals adopt the behavior above the discontinuous informa¬ 
tion transmission probability Xi. The horizontal white line is the 
critical degree exponent Vc, white circles and red dashed lines are 
simulated and theoretical results of Xi, respectively. Other parame¬ 
ters are 7 = 0.1 and po = 0 . 1 . 


nomenon in which the dependence of R{oo) on A can change 
from being continuous to being discontinuous. More partic¬ 
ularly, there is a critical degree exponent Vc below which the 
dependence is continuous [see region I in Figs.|3c) and (d)], 
while above the dependence is discontinuous [see region II 
in Figs. |3c) and (d)]. The value of can be gotten by bi¬ 
furcation analysis of Eq. ([T3 l. In region II, we also find that 
the discontinuous information transmission probability A^ in¬ 
creases with zz, since the fraction of hubs decreases with zz. 
The theoretical predictions of A^ can be gotten by bifurca¬ 
tion analysis of Eq. GSll, and the simulation results of A^ are 
predicted by NOI (number of iterations) method B. Re¬ 
gardless of network heterogeneity, our theoretical predictions 
about the behaviors of i?(c») have a good agreement with nu¬ 
merical calculations. The average relative error li^ between 
the two predictions of R{qo) for all the values of A and zz is 
less than 1 . 8 %. 

Einally, we study the effects of network topology on the 
final behavior adoption size R{oo) in Eig. |3for k = 2 and 
c = 1. We find that increasing zz can promote (suppress) 
behavior adoption at large (small) value of A. This phe¬ 
nomenon can be qualitatively understood in the in following 
ways Eor strong heterogeneous networks, the more 

hubs and a large number of individuals with small degrees are 
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FIG. 6 . (Color online) Behavior spreading on scale-free networks. 

(a) The final behavior adoption size R{oo) versus information trans¬ 
mission probability A for different degree exponents u = 2.1 (black 
circles), v = 3.0 (red squares) and u = 4.0 (blue up triangles). 

(b) The final behavior adoption size R{oo) versus A = 0.3 (black 
circles), A = 0.5 (red squares) and A = 0.8 (blue up triangles), 
respectively. The lines are the theoretical predictions from Eqs. © 
and O-Gll. Other parameters are c = 1, 7 = 0.1, k = 2, and 
po ~ 0 . 1 , respectively. 


coexisted. Since those hubs have more chances to expose the 
information, they adopt the behavior more easily even when 
A is small. However, the situations for individuals with small 
degree are just opposite. That is, the large number of individ¬ 
uals with small degrees hinder the behavior adoption for large 
value of A, thus, cause a smaller value of i?(oo). Through bi¬ 
furcation analysis of Eq. CSll, the dependence of R{oo) on A 


is always continuous for different v. 

V. DISCUSSION 

To study social contagion dynamics in human populations 
is an extremely challenging problem with broad implications 
and interest. For social contagions on networks, some inelas¬ 
tic resources (e.g., time, funds, and energy) restrict individuals 
to dedicate to social interaction, which have always been ne¬ 
glected in previous studies. In this paper, we first proposed 
a non-Markovian behavior spreading model with limited con¬ 
tact capacity, in which each adopted individual transmits the 
information to a fraction of his/her neighbors. We then de¬ 
veloped a heterogeneous edge-based compartmental theory to 
describe this model. The average relative error between the 
theoretical predictions and numerical calculations is less than 
1.8%. Through theory and simulations, we found that increas¬ 
ing the contact capacity c promotes the final behavior adoption 
size R{oo). With the help of bifurcation theory, we found a 
crossover phenomenon in which the dependence of R{oo) on 
A can change from being continuous to being discontinuous. 
More specifically, we uncovered a critical degree exponent i/c 
above which the crossover phenomenon can be induced by en¬ 
larging c. However, i?(oo) always grows continuously for any 
value of c when degree exponent is below i/c- 

Here we developed an accurate theoretical framework for 
non-Markovian social contagion model with limited contact 
capacity, which could be applied to other analogous dynami¬ 
cal processes (e.g., information diffusion and cascading fail¬ 
ure). Further more, how to design an effective strategy to con¬ 
trol the behavior spreading with limited contact capacity is an 
interesting research topic. 
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